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Abstract 

We propose a neutrino mass matrix model in which five neutrino species remain hght 
through the seesaw mechanism within a supersymmetric 3z/j^ + 3z//j framework. We con- 
struct such a model based on the nonrenormalizable terms in the superpotential con- 
strained by the discrete symmetry which may be expected in the models at the high 
energy scale such as superstring. We study the possible oscillation phenomena by fixing 
mass parameters so as to explain the solar and atmospheric neutrino deficits and also 
include a candidate of the suitable dark matter. We also discuss the charged lepton mass 
matrix based on this neutrino model. LSND results may be consistently explained within 
this model. 
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1 Introduction 



Neutrino mass is one of unsolved problems in the present particle physics, although it 
is a very important issue which can be a clue to go beyond the standard model. Ex- 
perimentally, in both particle physics and astrophysics there are a lot of indications for 
massive neutrinos |l|]. Analyses of observations of the solar neutrino and the atmo- 
spheric neutrino suggest the existence of neutrino oscillations. Especially, a recent Su- 
perKamiokande observation of the Zenith angle dependence of the atmospheric z/^ strongly 
suggests that i/^ oscillates to |Q . An interesting feature in these indications is that they 
may require wide range mixing angles, in particular, a maximal mixing among different 
neutrinos in addition to hierarchically small mass eigenvalues or closely degenerating mass 
eigenvalues. Although this smallness of masses may be considered to be explained by the 
seesaw mechanism in general, the hierarchy of masses and the mixing structure will be 
completely dependent on models. 

On the other hand, it has also been suggested that neutrinos with suitable masses could 
be a dark matter candidate for the explanation of astrophysical observations related to 
the structure formation of the universe. As such examples, we may list an active neutrino 
with (9(10) eV mass in the hot dark matter scenario (HDM) a sterile neutrino with 
(9(10-10^) eV mass in the warm dark matter scenario (WDM) and neutrinos with 0{1) 
eV mass in the cold + hot dark matter scenario (CHDM) |P, |^. When we consider the 
neutrino models, we may need accomodate this feature to the models. If we impose such 
a requirement, we can restrict the models in a rather severe way. 

Now various experiments using reactors, accelerators and underground facilities are 
proceeding and planned. In the near future these experimental results will be presented 
to inform us details of the neutrino sector and then the predictions for various neutrino 
phenomena including oscillations on the basis of possible neutrino models will also be very 
useful. Under this situation it seems to be a worthy and interesting subject to consider 
various types of model which can explain both of these neutrino deficits consistently and 
present a dark matter candidate. 

The introduction of a sterile neutrino is one way to this direction and a lot of works 
have been done by now llO|, ini, |12|, m. It has also been suggested that sterile neutrinos 



may play important roles in various phenomena [0, |T5|, 0. Although the existence of 
sterile neutrinos is an interesting possibihty, it seems to be rather difficult to find its 
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natural candidate in particle physics models. [| One of such reasons is that it is not so easy 
to introduce a substantial mixing between active neutrinos and light sterile neutrinos in 
the natural way [IB . 



The aim of this paper is to propose such a candidate and analyze neutrino oscillation 
phenomena on the basis of it. In this study the neutrino sector is extended into Su^ + 3z//j 
species, among which only five neutrinos remain light through the seesaw mechanism 
based on the mixings with a heavy right-handed neutrino [Np/). Two light right-handed 
neutrinos behave as sterile neutrinos. This seems to be natural from a viewpoint of the 
generation structure of other quarks and leptons, although the seesaw mechanism works 
in a different way from the one in the ordinary grand unified models such as SO (10). 
To realize the substantial mixing between these sterile neutrinos and active neutrinos we 
will consider the nonrenormalizable interactions coming from the fundamental theory in 
a high energy region. This may be considered as an example of the scheme proposed in 



Ref. |T8[. Phenomenologically, as partly discussed in Refs. [|I2|, |T3[, this model can give 
a framework to explain both deficits of the solar and atmospheric neutrinos consistently 
and also to present a candidate of the dark matter with a suitable mass for a certain kind 
of dark matter scenario. 

This paper is organized as follows. In section 2 we give a brief review on the oscilla- 
tion parameters to fix our notation. In section 3 we introduce our model and discuss its 
theoretical background. A possibility to realize this model in the super symmetric frame- 
work in the basis of nonrenormalizable terms in the superpotential is studied. Its detailed 
phenomenological analyses are presented in section 4. We show two typical parameter 
settings which have different features. Other oscillation processes than the ones related to 
the solar and atmospheric neutrino deficits are also discussed in each case. LSND results 
reported in Ref. []T^| may be simultaneously explained in this model. Section 5 is devoted 
to the summary. 



2 Oscillation parameters 

At first we briefiy review a basic formula for the neutrino oscillation. We define a mixing 
matrix V as Uf = J^a^af^a where (y'^)fa = V*j- and z/^ (a = 1,2,3, ■■■) is a mass 
^ For example, a candidate other than the right-handed neutrino has been proposed in Ref. p7|. 
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eigenstate. Uf (/ = e, /i, r, ■ ■ ■) is a weak eigenstate chosen in a way that both leptonic 
charged currents and a charged lepton mass matrix are diagonal. Thus V can be written 
as = [/('^)?7(')t ]-,y using diagonahzation matrices f/'-'-' and U^'^^ of the charged lepton 
and neutrino mass matrices M(') and M^^) 

In the following discussion, V is assumed to be real, for simplicity. Then V satisfies an 
orthogonality condition^: 

^»fVi3f = 6ap, J2 ^c.fVf3f = 6a(3. (2) 

/ 

Using mass eigenstates, a time evolution equation of neutrinos in a vacuum is given 

as 

where E is energy of neutrinos and stands for an a-th neutrino mass eigenvalue. This 
equation can be easily solved as 

i>„(t)=expUz^tjz>,(0). (4) 

By transforming this into a solution in terms of the weak eigenstates, the transition 
probability for Uf Uf during the time interval t is expressed as 



= E y'f'Vaf + 2 E E VafV^fV^rVf^f COS —^t , (5) 



where AM^^ = — M||. If we use a following relation derived from the orthogonality 
of V: 

E y^'^'f = -2 E E v^f,v^fVpf,Vpf, (6) 

" " /3(>q) 

we can obtain 



P.,^.^,(t) = -4^ Y: K/'K/V/^/'V^/^/sin^f ,p 

" /3(>a) V 



^"'^t\ . (7) 



^ We neglect CP phase here. In a model studied in this paper there are three charged lepton flavors 
and five light neutrino flavors. When we define V , it is necessary to extend t/^'-* into a 5 x 5 matrix by 
adding formally I's as diagonal elements. 
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Here we note that a contribution from a ct/^-sector to the i/f Vfi oscillation is 
represented by parameters (AM^^, — ^VafVafV^fV^f). In a sector where 14/ and Vg// 
correspond to the diagonal elements the amplitude ( or the mixing factor ) reduces to 
—AVaf'VfSf, as long as the mixings with other neutrino flavors are sufficiently small and 
Vaf ~ Vg/' ~ 1 is satisfied. In such a case these parameters can be understood as the 
usual two flavor oscillation parameters (Am^, sin^ 29). On the other hand, if the mixings 
with other flavors are not so small in the sectors specified by the off-diagonal elements 
Vaf and Vg//, there appear new contributions which are induced due to the existence of 
many neutrino flavors. They may be understood in such a way that an additional mixing 
between i/f^ and i/f^ induces the i/f <-> oscillation through the flavor mixings i/f — Vf^ 
and Uf/ — i/f^. Thus these processes are generally recognized as higher order effects in 
comparison with the direct two flavor oscillation concerning the mixing factors. Anyway 
it is necessary to be careful when we apply Eq. (7) to the results on the oscillation 
parameters obtained by the two flavor analysis of the solar and atmospheric neutrino 
problems. In particular, matter effects on the oscillations will not be observed in a simple 
way through an analytical study in the cases with many neutrino flavors. 

3 A model with sterile and active neutrino mixings 
3.1 Basic framework 

We consider a model containing three generation left-handed neutrinos Uf^ (/^ = e, /i, r), 
and right-handed neutrinos Nf^^ {fji — A,B,C). As a guiding principle to construct 
a neutrino mass matrix in this kind of models phenomenologically and reduce the free 
parameters systematically, we take a viewpoint that the realization of the mixing angle 
required in the solar and atmospheric neutrino problems is the most important clue. 
Along this line we first prepare mass terms which can produce such a mixing structure 
among the five light neutrinos through the seesaw mechanism. After this procedure we 
introduce mass correction to resolve the mass degeneracy in a consistent way with this 
mixing structure. 

Following this strategy, we require that six neutrino species have mass terms which 
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are written as 



5 5 _ ]^ _ _ 

-/:„,ass= E rrif f,^ fij f, + Y,mfijfNc + -MNcNc + h.c., (8) 
/,/'=! /=i 

where represents Uf^ and iVj^ (//? = A,B). Although the state identification is still 
not done at this stage, hierarchies among the above mass parameters are assumed asQ 



niff/ -C mi -C 777.2 -C ma ~ m4 -C m5 -C M. 



(9) 



Since mjj/ is small enough compared with others, we can neglect it for a while. As a 
result of the seesaw mechanism, a heavy right-handed neutrino Nc decouples from other 
neutrinos and a mass matrix for five light states ijjf becomes^ ||12| 



Mo = M 



2 

f^l A*l/^2 /^l/^3 /^lA*4 /^l/^5 
/il/i2 /ii /i2yU3 /i2Ai4 Ai2/i5 



yUl/^3 /^2yU3 /^3At4 yU3/^5 

2 

/il/i4 /i2/^4 /^3/^4 /^4Ai5 

,2 



(10) 



y yUl/i5 /^2/i5 /^3/^5 /^4/^5 / 

where fif = ruf/M. As is easily checked, Aio is a matrix with a rank one and diagonalized 
as U'^'^^ AioU^"^'^ by using the matrix 



O 
1 





/i2 


/ii 












6 


6 






/Xi/i3 


At2/i3 









66 


6 






/il/i4 


/i2/^4 


/i3/i4 


6 




66 


66 


66 


6 








/W3/^5 


/i4/i5 




66 


66 


66 


66 






/^2 


/^3 


/i4 


V 


6 


6 


6 


6 










6 
6 / 



n+l 



where ~ X] /^/' is an undetermined 4x4 matrix which should be introduced because 

of the mass degeneracy. 

The inequality between mi and TO2 is reversed from the one in Ref. [ p^ . The present one should 

be used to make the MSW mechanism applicable in the {ipi, ^2) when the state identification is assumed 

such as T/^i = and ^2 = i^e which will be used in this paper. Some related equations in Ref. 

should be replaced by the ones presented in this paper. 

It should be noted that we approximately call these states as ■)/'/• They are not pure V^/'s but have 

a small mixture component from Nc- 
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To resolve this mass degeneracy and fix U^'^'' completely, we now switch on rather 
complicated mass corrections rriffr to the five light neutrinos in order to yield the hierar- 
chical mass eigenvalues without disturbing the mixing structure U^'^^ . As such masses we 
consider the simplest example as 





' Alii 


B^l^2 






JF/ilyUs ^ 






Cfil 
























Sfl2fi4. 




Vfil 


T /i4/i5 














JF are parameters which 


satisfy 








A- 


V B- 


V 








B- 


V C- 


V 


/if 


V + £ 


i4' 



(12) 



(13) 



This matrix can also be checked to be diagonalized by t/^'^-' and the mass eigenvalues of 
A^o + -^per is obtained as 

Mi = {A- V)fil, M-2 = 0, M-s = {V- S)fil Mi = {V + S)fil Mg = Mfij. (14) 

Hereafter we use numerical indicies with a bar to specify the mass eigenstates. A 55- 
element of A^per can take its value in a range like \Q\ ^ 0(|JF|) to guarantee the approx- 
imate diagonalization by U^"^ because the fifth mass eigenvalue Mg coming from A^o is 
rather large. After the introduction of this correction the mass degeneracy is completely 
resolved and O in Eq. (11) is determined as O = 1. 

Here it is useful to note how many number of parameters are included in this model. 
As easily found from Eqs. (11) and (14), there are nine free parameters. They may be 
taken as /i/(/ = 1 ~ 5), M,C,V and S. It is convenient for the later study to present 
concrete expressions of oscillation parameters in Eq. (7), in particular, the amplitude 
—4:Vaf'VafV^f>V^f by using these parameter^. For the ipf ipf oscillation their expres- 
sions due to the a/5-sector are presented in Table 1. Although there are small negative 
contributions to V^^,^^^ from some a/5-sectors, we omit them from Table 1. The detailed 
explanation of numerical values listed in the columns (I) and (II) of Table 1 is given in 
the next section. 



For a while, we confine our attention to assuming that a charged lepton mass matrix is diagonal 
so that U^^^ = 1. In section 4.3 we extend our analysis to more general cases. 
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3.2 Construction based on nonrenormalizable terms 

In the previous subsection we introduced our mass matrix A^o + -^per in the phenomeno- 
logical way so as to reahze the substantial mixing between sterile and active neutrinos 
with nondegenerating masses. However, we can also construct it under a certain theoret- 
ical background. We present such an example in this subsection. To make our argument 
definite we consider a certain type of supersjTiimetric nonanomalous extra U(l) models 



with a D-flat direction for the extra U(l), which comes from Eq unification group ||20|| . 
This kind of models may be expected to be induced from the perturbative superstring. 
The following argument will be straightforwardly applied to other type of models, for 
example, anomalous U(l) models. 

This example contains three right-handed neutrinos Nf^ and pairs of singlets {Se, Se) 
{£ = 1, 2), where Nf^ and Se belong to 27 chiral superfield of Eq and a conjugate partner 
Se belongs to 27* chiral superfield. The extension to the case with £ > 3 is straightforward. 
Nfjf^, Se and 5*^ are singlets under the gauge groups of the standard model but have the 
extra U(l) charge. Nf^ and Se are assumed to have the same charge of this extra U(l). 
The gauge invariant relevant terms in the superpotential to the present investigation are 
composed by two parts, 

Wo = E ^^siSeSef'^ (15) 

£=1,2 ^^^G 

f rJ' 

fr.fn Mn 



IlJr ^VIg 

+ E {S,S,y^'^- (5,^,)"^^^^-^^^^^ SlH^,L,,Lr^, (16) 



IlJ'^ '^'^ 

where the parameters determining the power structure of each term are integers and 
satisfy 

Ke>2, Pf,r^-2>pf^r^>0, Qf,f,>qf,f,>0, i?/,/^ > r^,;, > 0. 

Mg is a suitable scale such as a string scale or a Planck scale. The higher order terms can 
be neglected in comparison with these lowest dimensional terms. The scalar potential for 
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the singlet scalars S(, and Sg, can be written as 

y = 'fi:Qi{\Si\'-\sA'f 

^ e=i,2 

The first hne is a D-tcrm contribution of the extra U(l) and the second hne represents 
an F-term contribution coming from Wq and soft supersymmetry breaking scalar masses 
which we assumed as = = m^, for simplicity. Clearly this scalar potential has 
D-flat directions l^*^! = |^^| = U(, along which the potential minimum is realized. After 
soft supersymmetry breaking masses are introduced as in Eq. (17), the magnitude of the 
intermidiate scale is determined by the nonrenormalizable F-term contribution in the 
second line as 



/I \ 2X^-2 

where the value of ug, is crucially dependent on an integer Kg^. Moreover, as easily seen 
in Eq. (16), depending on the values of powers in the nonrenormalizable terms such as 
Pfnf'^, QfLf'a ^^"^ ^fLfL-i various types of neutrino mass terms can be induced through the 
effective couplings controled by the power of uiIMq after the vacuum expectation value 
ui becomes nonzero. 

Since there is no restriction on the right-handed neutrino mass terms due to the gauge 
symmetry unlike the usual Grand Unified models, this framework can contain light sterile 
neutrinos in addtion to the ordinary three active light neutrinos. If we take the basis -0/ 
in Eq. (8) as t/jf — {Na{= i^sj, i^/^, Nb{= i^s^)):^ the mass parameters in Eq. (8) can be 
written by using the parameters in Eq. (16) as 

where we use the definitions e = Ui/Mq and 5 = U2/U1. In Eq. (19) the indicies 
/l = e,fj,,T should be interpreted to correspond to the numerical indices 2, 3 and 4 
in Eq. (8). Then the parameters /// are also expressed as 

I.AC UlC 
- 1S_ 2iPAC-Pcc)fj2{pAC-pcc) - 1d_^ 2iQ^^c-Pcc)+3x2iqfj^c-PCc) 

^1 ~ J^CC ^ " ^ f^fL — f^CC u ' 

s s ^ 



^The phenomenological validity of this identification will be discussed in the next section. 



9 



^5 = %_,2(P,c-Pcc)52(p,e-Pcc). (20) 
kg 

Moreover, based on the superpotential Wi we can write the neutrino mass matrix A^per 
under the same basis as follows, 



J^fLA^2Qf^A^2qf^A^^ k{}^e^^fLBS^1fLBy^ , (21) 

y kj^e^P^B-^d^P^^Ui kfj^^e^^iLB^^^fLBy^ J^BB^2PsB-SS^PBBu^ J 

where f 2 is a VEV of the doublet Higgs scalar Ii2- Major ana masses M. /^/^ can be caused 
by the last terms in Eq. (16). However, even if Rfj^f^ = 0, these Majorana masses are 
~ 10~^ eV and then generally too small for the explanation of the solar and atmospheric 
neutrino problem 0]. Within the interactions included in Eq. (16) M.f^f' can be also 



induced as the mixing terms composed by the first and the second terms of Eq. (16). 
But it is difficult to arrange in the way to satisfy the condition given by Eq. (13) and 
guarantee the sufficient largeness of eigenvalues. Finally we are forced to introduce new 
interaction terms. Here we introduce a triplet Higgs scalar $ and consider the following 
interaction terms in the superpotential [|: 

= ^$^(S2^2)*^-^i(5i5i)^^-^l"*^'-^l<l>L;,L;, . (22) 

This gives the Majorana masses 

Mf^f^ = k!r'^'^e^^f^''L6''^f'LVT, (23) 

where vt is a VEV of $. There is a constraint on vt from an electroweak p parameter. 
To satisfy its constraint it will be necessary to take it as f t ^ 1 GeV. 

What kind of the mass matrix A^per is induced in this model is completely dependent 
on the values of the lowest powers and couplings of each terms in the superpotential shown 
by Eqs. (15), (16) and (22). In order to constrain the superpotential, it is necessary to 
be able to introduce a certain kind of symmetry which can forbid the lower dimensional 
terms consistently. In many works an Abelian horizontal symmetry has been used 



to constrain the nonrenormalizable superpotential, which can induce the small neutrino 



^ If we try to construct the model including a triplet Higgs scalar in the perturbative superstring, it 
must be necessary to consider the models with higher Kac-Moody level. 
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masses with a favorable texture after the breakdown of this horizontal symmetry due 
to the VEV of some scalar fields. In our model the discrete symmetry may play the 
same role. To realize the superpotential starting from high dimensional terms, we need 
rather complicated higher order discrete symmetry. Here we assume Zg x as a concrete 
example of such a discrete symmetry and take the charge assignments for the relevant 
fields under this symmetry as the ones shown in Table 2. The charge assignment for this 
discrete symmetry allows the lowest terms with Ki = K2 = 9 in Wq. This results in 

~ ~ 10^^ GeV and then e ~ 0.1 and 5^1. We should note here that in this case 
e is not so small that the extremely high dimensional terms are necessary to realize the 
desirable neutrino masses. If we adopt a smaller and make e small enough, we do not 
need such high dimensional terms but additional fine tunings for coupling coefficients are 
required to induce a necessary mass pattern. In such a case the main feature of the mass 
pattern is determined by the tunings of the couplings. Thus we do not take this way here. 

Terms in Wi and W2 are also constrained by this discrete symmetry and the following 
values for the powers in Wi and W2 are allowed as the lowest ones: 
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The present examples may not be considered to be realistic because of their extremely 
large dimensions of the necessary nonrenormalizable terms. Although this feature seems 
to be general in this type of models, the situation may be made mild to some extent by 
considering a type of extra U(l) models in which a role of S^Si is replaced by suitable 
elementary singlet fields with D-fiat directions. Anyway this example shows the direction 
how we can constrain the nonrenormalizable terms in the superpotential. 

Next we study the neutrino mass matrix induced by this superpotential. If we take 
V2 ~ e^^Mi, vt ^ i GeV, these values give Hf and A4per as follows: 

jll — Kg e , II2 — Kd ^ ^ l^'i — I^D ^ ^ — t^D ^ : A*5 — % ^ , 
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UAA 

h.eA 

heChAC 



UAC 
''D 
UtA 

tLtC Lac 

Kjj Kg 
LAB 

rv Q 



\ h.BC uAC 
\ Kc Ave 



LeA 
LeCuAC 



Lee 

1^2 



{kff 



h.eB 



UeCuBC 
Kjj Kg 



It^ClAC 
Kjj Kg 

Leu 



LeC Ul^C 
Kq Iv£, 

nj'-p 9 

i4 
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hM<^ utC 



LiJ-CuBC 
Kjj Kg 



1^31^4 



UtA 



I.TCI.AC 
Kjj Kg 



her 
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Kjj Kjj 
Lt^T 



Uf^CurC 
UTT 
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UtCuBC 
Kjj Kg 



I.AB 

LbcLac 

Kg Kg 
UeB 
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LU-B 



I.BC 
''D 

h.TB 



h.TC h.BC 
tv^ Ivg 

Lbb 



/i4/i5 



(25) 



where k^-^' = k^-^' /kg^. To obtain the phenomenologically favorable hierarchy among 
Hf {f = 1 5) and reahze the mass matrix form defined by Eqs. (12) and (13), we need 
additional hierarchical structure in the coupling constants k^-^' ■ Among these tunings of 
coupling constants we introduce the following useful parameters expressing the freedom 
which can not be determined by the present experimental results in our state identification 
used later: 

« = 2 - log.o , ^'=-logio(^)- (26) 

Using these parameters and if we make the suitable assumptions on k^f , fj, f can have the 
following hierarchy: 



^12+a 
/^2 ~ e , 



^12 

AI3 ~ e , 



//5 ~ e 
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(27) 



Although there arc many possibilities for A^per depending on the assumption on the 
coupling constants, it will be useful to present examples to see what kind of tunings of 
the coupling constants are required to construct the interesting A^per- Here we give two 
examples: 



(I): 
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(II): 



e //3//4 



-e//i/i5 



(28) 



e/ii/i4 e/i2At4 eAi3yU4 

\^ -6/^11^5 -61^21^5 -(^1^31^5 -6/^4/^5 

These examples have the interesting penomenological features as explained in the next 
section. Although we need additional tunings for the coupling constants to satisfy the 
condition shown by Eq. (13) as found from these examples, such tunings seem not to be 
so hard but rather mild as we can find by comparing Eqs. (25) and (28). We believe that 
it does not spoil the interesting feature of the present model. 

For the charged lepton sector the same discrete symmetry can also determine the 
structure of the nonrenormalizable terms in the superpotential 



E 



k 



^2Xf 



-{S2S2Tf^f^{SiSiffLfR-^fLfRLf^HiEf^ 



(29) 



fTTfn Ma'-'- 

following the charge assignment in Table 2. By the use of this symmetry we can write 
down the charged lepton mass matrix M^^^ as follows: 



^ kl^vie^ kfvie^ kfvie'^ \ ( cl ^ 



{gr, Hr, Tr) 



(30) 



k'^v^e^ k'^v-^e^ k'^Vie'^ 
\ k]!vie^ kl^v^e'' k][v,e^ ) 

where coupling constants again need to be tuned to realize the desirable charged lepton 
mass eigenvalues. Here we assume that k^ < k]^ and also the off-diagonal couplings are 
smaller than these diagonal couplings. Then this mass matrix can be diagonalized by the 
bi-unitary transformation and results in the eigenvalues 



rrie ~ k'^Vie , ~ k^ Vie , m,- ~ k'^Vie . 

The diagonalization matrix C/^'^ is almost diagonal and can be written as 



(31) 



1 -{kt/ei)^ -{k%^iki-y 



(32) 



^ {k'^/k'^)e^ k'^ /k][ 1 j 

A possibility of the mixings among sterile neutrinos and active neutrinos has already 
been proposed based on the nonrenormalizable interactions in the superpotential in the 
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context of superstring inspired models in Ref. [0. Our model may be considered as a 
concrete example of its realization. Although our scheme may need some tunings at least 
for the coupling constants in the superpotential to make the mass matrix the required 
form defined by Eqs. (12) and (13) in the precise way, it is interesting that only the above 
values of -P/jj/^, QfLfn Tj^/^ can approximately make it close to the required form. If 
we assume ui ^ U2, more available freedom may be applicable. This kind of models may 
be recognized as one of many candidates for the possible neutrino mass matrix realized 
in the promising supersymmetric models inspired by perturbative superstring. 



4 Analysis of various oscillations 

We apply our model to the analysis of neutrino oscillations. Deficiencies of the solar 
neutrinos and the atmospheric neutrinos ^ have been suggested to be explained by 
i^e ^ i^x and z/^ ^ Uy oscillations, respectively. Within a two flavor oscillation framework 
the neutrino squared mass differences and mixing angles predicted from these observations 
are for the solar neutrino problem^ p|, ^ , 

Am^^^^ ~ (0.3 - 1.2) X 10" sin^ 29 ~ (0.4 - 1.5) x 10"^ (33) 

and for the atmospheric neutrino problem 0, ^ , 

Aml^^^ ~ (4 - 6) X 10'2eV^ sin^ 29 > 0.85. (34) 

We also take account here that the existence of one neutrino species with mass such as 
1 ~ 10^ eV has an interesting relevance to the astrophysical observations for the large 
scale structure of the universe. 

In order to present the consistent explanation for the neutrino mixings shown in Eqs. 
(33) and (34), we need to identify the five light states tpf with physical neutrino states. 
In this consideration the constraint from the standard big bang nucleosynthesis (BBN) 



23[| may be useful because it can constrain a sterile neutrino sector. The BBN predicts 
that the effective neutrino species during the primordial nucleosynthesis should be less 
than 3.3. This fact severely constrains the mixing angle 9 and the squared mass difference 
Am^ between a sterile neutrino (u^) and left-handed active neutrinos which mix with it 



* There are also large mixing solutions. However, we do not consider these solutions in this paper. 
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p^ , As long as we do not assume the large lepton number asymmetry at the BBN 
epoch |j26|, these constraints rule out the large mixing MSW solution of the solar neutrino 
problem due to z/g Us and also the explanation of the atmospheric neutrino problem by 
^s- Taking account of these facts, we concentrate our study on a possibility such 
that ipi and are right-handed sterile neutrinos and Vg^, and il)2-, ips and ip4^ are active 
neutrinos z/g, z/^ and z/^. In this identification the solar and atmospheric neutrino deficits 
are considered to be explained by the small mixing MSW solution due to u^. — > z/^^ and 
the z/^ ^ Ur oscillation in the vacuum, respectively.^ In Table 1 this identification has 
been assumed. 

Now we impose that the two flavors oscillation scheme is good enough for these oscil- 
lation processes. And then we can determine some of the mixing parameters numerically 
based on both oscillations in the (z/s-^,z/e) sector with AMj^ and in the (z/^, z/^) sector 
through a mode with a = 3 and /3 = 4 which are shown in Table 1. As f/*^'-* = 1 is 
assumed here, the desired mixing angles in Eqs. (33) and (34) can be realized by setting 

16 < ^ < 32, 0.44 < ^ < 2.3. (35) 

In the following discussion we take these values as /i2//^i — 25 and fis/fi^ — 1, for sim- 
plicity. To investigate other oscillation processes it is convenient to use the parameters 
a and b introduced in Eq. (26) which gives their physical meanings as the ratio of cou- 
pling constants. They can be also expressed as /i2//i3 = 10^" and n^l = 10^^ and 
parametrize the mixing among different neutrino species. There seem to be no quantita- 
tive constraints on these parameters at the present stage. From the viewpoint that the 
solar and atmospheric neutrino deficits are explained by the two flavors scheme for the 
oscillation processes, it is enough for them to be sufficiently large. 

However, if we take account of the BBN constraint in the more quantitative way, the 
a and h dependence of the mixing parameters seems to allow us to restrict the region of a 
and h. As mentioned earlier, the restriction on the number of the effective neutrino species 
during the primordial nucleosynthesis gives the condition on the i's2~^e,fj.,T sector. Here it 
is sufficient to consider the most stringent one. In the two flavors oscillation scheme when 
^ There is another possibihty that the solar and atmospheric neutrino deficits are explained by the 
f/j and Vfj, ^ Vt oscillation, respectively. However, in such a case one light sterile neutrino plays no 



role in the oscillation phenomena and it can be reduced to the model considered in Ref. |13 
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m,, > m,, is satisfied, it can be formulated as 

Am^ sin^ 2^ < 3 X 10~^ eV^ for iu^,,r,i^s)- (36) 

If we can apply this constraint to the (i/^, i^sj) and {ur, 1^32) sectors in Table 1, nontrivial 
constraints on b can be obtained.^'' This condition may be written as 

10-^\Ml - Mi) < 10-^-^ (37) 

If we require that the heavier right-handed neutrino z/^j can be a dark matter, a value of 
Mg should be fixed so as to be a suitable value and then the parameter b should satisfy 
the following conditions: 

b > 1.6 (Mg ~ 1 eV), b > 2.1 (Mg ~ 10 eV), b > 2.6 (Mg ~ 10^ eV). (38) 

If we consider the {I'si, ^e,^i,T) sector with m^^ < ml^^^, through the BBN constraint 
given in Refs. |]2^, ^ there appears no condition on a for the case of AM^ ~ 10"^'^ but 
a ^ 1 should be satisfied for the case of AM^ ~ 1. The situation completely depends on 
the suqared mass difference among and i^e,^i,T- Of course, if we assume the existence 
of the large lepton asymmetry at the BBN epoch, these constraints can disappear p6|l . 



However, a lower bound on the parameter a can be always obtained from the condition 
on the amplitude in Eq. (7) such as —^VafVafVpfiVpf < 1. Using Table 1, a weak 
constraint on a is brought as a 0.15 from this condition. 

In order to realize the desired squared mass differences in Eqs. (33) and (34), we need 
to require 

Ml ~ 10"^-^ eV, |M| - M|| ~ IQ-^-^ eV^. (39) 

If we take Mg ~ 10 eV as an example and use the first one in Eq. (39), the parameters 
in our model defined by Eq. (8) should be settled by using Eqs. (14) and (39) as, 

mi ~ 10-^-^~('^+'')M^ ma ~ lO^^^^'^+^'^M^ 

We should be careful in this application since this constraint has been derived in the two flavors 

oscillation scheme. The BBN constraints are crucially affected by the interaction with the plasma at 

finite temperature so that the situation may be changed in many flavors case from the one of two flavor 

oscillation scheme. Athough we need a numerical calculation for the correct analysis of this aspect, such 

a study is beyond the scope of this paper. However, this kind of consideration may be useful to show the 

importance of the BBN constraint in the model building. 

^^This value should be changed dependeing on what kind of dark matter scenario is considered. Here 

we take a conservative value not far from required values in the various models |^, 0, ||] • 
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m^^m^^ IQ-^-^Ma, mg ~ 10" (40) 

where we use a GeV unit, and a and h should satisfy the constraint given by Eq. (38) 
and a ^ 0.15. To make |M| — M|| a value presented in Eq. (39), there can be many 
possibilities for the values of Mg and Mj. Here we consider the following typical two 
cases : (I) Mg < ~ lO'^-^ eV, and (II) ~ (> lO'^-^ eV). In the case (II) we 
take two eigenvalues as Mg ~ Mj ~ 1 eV which has been studied in the various works 



| T0[] as an interesting example, although such a choice requires a rather strict degeneracy 
between the third and fourth mass eigenvalues. Numerical expressions of the oscillation 
parameters for each of these two cases are given in the columns (I) and (II) of Table 1. 
We should also note that these cases with certain values of a and b can be realized as the 
two models (I) and (II) presented in the previous section. 

If we observe Table 1 taking account of Eq. (38) and a ^ 0.15, we immediately find 
that the very restricted oscillation modes can effectively occur and others are negligible 
because of the small amplitudes ( mixing angles ). There are two a-independent processes 

(z/,,,z/e) with AMl^, (i/^,z/,) with AMfj, (41) 

and also as the a-dependent but non-negligible interesting oscillation modes, we have 

{u,,iy^) with AM|3,AM|4, (z/^, z^.) with AMf^. (42) 

As already mentioned, two processes given in Eq. (41) can be treated within the two 
flavors oscillation scheme, as long as the parameter a takes a suitable value which can 
guarantee such a treatment. In that case they can be used for the explanation of the 
solar and atmospheric neutrino problems in both cases of (I) and (II). On the other hand, 
although some of the processes listed in Eq. (42) come out as the effects due to many 
flavors existence, they may bring about the important contributions to the tpf tpf/ 
oscillation according to the value of a. Next we examine both cases (I) and (II) in more 
detail and also discuss these processes in each case. 

4.1 Case I : M 5 < - IQ-^ ^ eV 

In this case [0, from Eqs. (13), (14), (39) and (40), we should take the parameters as 

A ~ io2("+^)-o-7m, B ~ -102("+^)-3-5m, 

C ~ P ~ £ ~ 10^''"^-^M, IQ-'^-^M, (43) 
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where M5 ~ 10 eV and a > 0.5 are assumed. The allowed range of a can determine the 
sign of B, since B < requires {B — T>)/T> < — 1 and then a > 0.5. V ^ £ should be 
reminded in order to resolve the mass degeneracy = M3 = 0. The consistency of 
this model requires that the largest mass parameter A should be smaller than Mpi. This 
requirement brings about an additional constraint on the parameters a and 6, 

a + 6<0.35 + logio(^)'. (44) 

If we take M ~ 10^^ GeV which is realized in the examples defined by Eq. (24) and use 
Eq. (38), we can constrain the value of a, for example as follows: 

a < 1.75 for h > 2.1 (Mg ~ 10 eV). (45) 

It should be also noted that we can judge the phenomenological validity and consistency 
of the constructed models based on the discrete symmetry such as the ones shown in Eq. 
(28) by comparing it with Eq. (43) and studying whether a and h satisfy the conditions 
given by Eqs. (37) and (44). In the model defined by (I) of Eq. (28) we obtain a ~ 1.5 
and 6 ~ 2.2, which satisfy these conditions. 

In the basis of these knowledge we can readily investigate the processes shown in Eq. 
(42). As easily found in Table 1, if we take a ~ 1.25 and Mg ~ 10"^'^, the oscillation 
parameters of the first process with AMlg in Eq. (42) also seems to take appropriate 
values for the small mixing MSW solution of the solar neutrino problem. In this case 
the light sterile neutrino may not be necessary for the explanation of the solar neutrino 
deficit. For smaller a, anyway, this process with larger Mg, (i/g, i/^) with AM|j and (i/g, v^.) 
with AMgj may be good targets of long basehne experiments. However, these processes 
in the interesting regions of Mg and a does not seem to be the ones suitable for the 
two flavors treatment. Especially, in these processes the matter effect seems not to be 
analytically estimated in the precise way and we may not be able to apply the condition 
given in Eq. (33) to this case naively. The numerical analysis for the oscillations among 
many flavors will be indispensable for more detailed study. If we want to guarantee the 
validity of the two flavor oscillation analysis for both the solar neutrino and atmospheric 
neutrino problems due to <-> Ve and <-> we need to require a ^ 1.3 and then 
M < 10^^ GeV, which can be satisfied in our model defined by Eq. (24). In such a case 
all other oscillation processes than and i/^ ^ listed in Eq. (41) unfortunately 

seem to be inaccessible experimentally. 
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In this case all active neutrinos are too light to be a hot dark matter but z/^a may be 
a warm dark matter with m^^^ = 0(10 — 10^) eV and Q,^^^ ~ 1 0. The main problem is 
how it is sufficiently produced at the early universe. There seem to be two possibilities 
for its production. If Ug,^ has an interaction with other light fields to be in the thermal 



equilibrium and then decouples as a relativistic particles, there is a relic 



n..,h^ = -^, (46) 

2 5f* 8.5 eV 

where h = iJo/(100km/sec/Mpc) and Hq is the present Hubble constant, is the effective 
degrees of freedom of the light fields at the decoupling time. In the present model i/^j 
has the interaction with other fields due to the extra U(l) gauge symmetry which breaks 
down at a very high energy scale ui. The above formula is applicable to Us^- If we assume 
h ~ 0.5 and g^, ^ 300 which is usually expected for this type of supersymmetric models at 
the decoupling epoch in the very high energy scale, we have ~ 1 for Mg ~ 10^ eV. If 
inflation occurs after the decoupling of u^^, however, this possibility cannot be applyed. 

Another possibility is the production through the z/^,r-'^s2 oscillation as suggested in 
Ref. 0. In this case z/^,r are at ffist in the thermal equilibrium and then decouples 
satisfying the same relation between and m,^ given by Eq. (46). During this period 
is considered to be produced from z/^,- through the oscillation process u^^r ^ ^S2- When 
we take this possibility, the constraint on Mg comes from the requirement for both of 



the sufficient abundance and the consistency with the BBN |^5[. If we follow Ref. 0, 
the ratio of the distribution functions fg and Ja of sterile and active neutrinos can be 
estimated in our model as 

f^^Q^ (m^\^ / keV \ ^ Mg^ 

/a gl'^\ eV ) [Mf^lJ' nA M-Ja ^ ^ 

By applying the hot relic relation given by Eq. (46) to these formulas and remembering 
Mg = M/ig, we can derive the following relation: 

Mg ^ = 0.22/inif (4g) 

On the other hand, the BBN constraint requires jsj ] a ^ 0.4 and if we impose this on the 
latter of Eq. (47) and use Eq. (46), we can obtain 

Mg > 230/i2r]^^ eV. (49) 
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In this case we should take ~ 10.8 and then the required values for Mg and h = 
In ifi^/fis) are Mg > 58 eV and b > 2.7 for Q^^^ ~ 1 (WDM) and Mg > 17 eV and b > 2.4 
for Q^^,^ ~ 0.3 (CHDM) |g. This value of Mg for the case of CDHM is so large that 
the free streaming length of becomes too short and then seems not to contribute the 
structure formation at the supercluster scale 0. Thus z/^j can be expected only to play 
the role as the warm dark matter. In that case a cannot take large value and a ^ 1.1. 

4.2 Case II : Mg ~ M4 ~ 1 eV 

In this case we need a rather strict fine tuning like Mj — M3 ~ 10"^''' eV. If we assume 
such a fine tuning, the parameters can be settled as 

A ~ io2(»+^)-o-7m, B ~ -102("+^)-3-5m, 

C ~ P ~ lO^^-^M, S^IO'^^-^M, J^r^-IQ-^M (50) 

where Mg ~ 10 eV and a > 1.3 is assumed. The condition given by Eq. (44) should be 
satisfied also in this case. If we try to realize this case in terms of the model defined by 
(II) of Eq. (28), we obtain a ~ 1.5 and b ~ 2.2 and the required conditions for a and 
b are fulfilled. The difference from the case (I) concerning the oscillation phenomena is 
that the rather large squared mass difference such as 0(1) eV^ can appear between Uf, 
and Ufj^^r- 

For such a squared mass difference, we can expect that the processes in Eq. (42) 
become interesting ones from the experimental viewpoint. If we apply the results of BNL 



E776 p8| and KARMEN p9| experiments for the Ue appearence through z/^j to the 
first process in Eq. (42), we can obtain a new lower bound on a. Two fiavors oscillation 
analysis of the data obtained by these puts the most stringent bound on a mixing angle 
among i/g and such as sin^ 2^ ^ 7 x 10~^ for Am^ = 0(1) eV^. Using Table 1, we 
can obtain a ^ 1.23 from this bound. This also satisfies the constraint from the z/g 
transition from Burgey [^] which requires a ^ 1.0 for the first one in Eq. (42) and also 
the constraint from the BEN which was mentioned before. 

In this context the interesting experimental results are the ones of LSND [TU|. The 



evidences for the oscillations and —* Ve have been reported there. One of nice 

features in this case is that these LSND results seem to be explained in the present model 
in terms of the first two processes in Eq. (42). In fact, as we take AMI3, AM|j ~ 1 eV^, 
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for this squared mass difference the LSND results require that the mixing angle should be 
3x10-3 < sm'^29< 1.5 x lO'^ for z/^ ^ z/^ and 1.5 x lO'^ < sin^ 2^ < 1.5 x IQ-^ for z/^ ^ z/g. 
These can be realized by taking 1.1 ^ a ^ 1.4 and 0.6 ^ a ^ 1.6, respectively. If these are 
combined with the above results of BNL E776 and KARMEN, a can be constrained to 
a very narrow region 1.23 ^ a ^ 1.4. If a takes a small value which cannot explain the 
LSND results, the second process in Eq. (42) may become a very interesting target in 
the future experiments. An interesting aspect of this case is that the above region of a 
relevant to the LSND results can put the constraint on the value of Mg which is the mass 
of a candidate of the dark matter. In fact, if we assume a ~ 1.4 and M ~ 10^^ GeV as 
an example, we obtain the bound on Mg as Mg ^ 36 eV by combining Eq. (44) and the 
BBN constraint given by Eq. (37). 

Anyway this case in our framework corresponds to an interesting realization of the 
model which has been pointed out by various authors |]10[, in which the solar and at- 
mospheric neutrino deficits and the LSND results can be simultaneously explained and 
additionally and fir can be the hot dark matter candidates in the CHDM scenario. In 
this case we may naturally ask the physical role of i>s2- The first interest is whether i>s2 
can have some affection for the structure formation or not. We can estimate this by using 
Eqs. (48) and (49) as in the previous case. For example, if we assume Q^^^ ~ 0.1, we 
have Mg ^ 5.8 eV and also b ^ 2.2 which is included in the allowed region presented by 
Eq. (37). This suggests that may have some affection on the structure formation as 
a part of hot dark matter in the CHDM scenario. Another interesting possibility of the 
physical role of h's2 niay be an effect on the leptogenesis discussed in Ref. [|1^]. We do not 
study it here but it may be an interesting aspect of our model. 



4.3 Relation to the charged lepton sector 

It may be useful to comment on the constraint on the charged lepton mass matrix in the 
present framework. Although we have assumed that the charged lepton mass matrix is 
diagonal up to now, it is also related to the neutrino oscillation phenomena through the 
mixing matrix V as shown in Eq. (7). If we once fix the mixing matrix of the neutrino 
sector as Eq. (11), the neutrino oscillation data can constrain the structure of the charged 
lepton mass matrix. 

Here we consider two typical examples as the charged lepton mass matrix. The first 
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example (A) is represented by Eq. (30) which is obtained in the basis of the discrete 
symmetry. Using the mixing matrix given in Eq. (32), we can obtain the mixing 
matrix elements Vaj defined in section 2 as follows, 



1 u e*"^ e*^ 

Vie ^ , Vj. ~ Vir ^ (51) 

where we have introduced phases a and r for completeness. 

As an another typical phenomenological example (B), we adopt a Fritzsch mass matrix 



3111 for the charged lepton sector. Although it is not relevant to our construction of the 
neutrino sector based on the discrete symmetry given in Table 2, it will be useful to find 
the feature of the mixing matrix presented in Eq. (11). Using a well-known formula in 
the diagonalization of the Fritzsch mass matrix as we can get the mixing matrix 
elements for the lepton sector asQ 



V rn^ /i3 /i3 y 



V,„.^ e-e-M^ . V,.^ e-J^ + e-M. (52) 



where me,m^ and m,- are charged lepton mass eigenvalues. 

In both cases V is found to have the similar form except that the latter example has 
extra contributions to the off-diagonal part compared to the for mar one. We present 
numerical values of the representative oscillation parameters for suitable settings of a and 
b in Table 3. From this table we find that both solar and atmospheric neutrino deficits 
can be simultaneously explained for these charged lepton mass matrices. Related to this 
point it may be useful to note that these charged lepton mass matrices have no effect on 
Vie and Vgsj, which are the elements of the above mixing matrix relevant to the ^ 

These expressions are somehow different from the ones given in Ref. jl^ since we take the mass 
hierarchy given by Eq. (9) which should be assumed for the explanation of the solar neutrino problem. 
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oscillation as found from Eq. (7). This feature is very different from the models in which 
the similar mass matrices are assumed for both the charged lepton and neutrino sectors. 
In those models, the mixings Vg^ and Vgg has a tendency to become too large due to the 
contribution from the charged lepton sector without assuming suitable values of phases 
to explain the solar neutrino deficit by the small mixing MSW solution for if we 

keep Vg^ and to be suitable for the explanation of the atmospheric deficit due to 



l^r 



32| , |33| . This aspect also appears in the {Uf,, z/^) mixing of the case (B) in Table 



3. The present model does not suffer from this problem as a direct result that the solar 
neutrino deficit is explained by the z/g z/^^ oscillation due to the introduction of a sterile 
neutrino. These two examples show that as long as the U^''^ is approximately diagonal 
in the similar way to these examples, our scenario is always expected to be applicable 
independently of the details of the charged lepton mass matrix. 

When we adopt these charged lepton mass matrices, the LSND results can be also 
explained in the case (11) in the same way as discussed in the previous part as long as 
we take a in a suitable region around a ~ 1.3. The situation on the consistency with the 
BBN constraint is also similar to the case of [/W = 1 since the charged lepton sector has 
no important effect in the (z/gj,//^) sector. No contradiction happens against the BBN 
constraint through the oscillations in the (us^jI^t) sector if b and Mg are in the suitable 
region shown in Eq. (38). This feature can be seen in Table 3. For other processes the 
same discussions presented in the present section are also vahd in the present case. 



5 Summary 

We proposed the neutrino mass matrix in the Su^ + 3z/^ framework, which could be 
constructed using the nonrenormalizable terms in the superpotential constrained by the 
suitable discrete symmetry. We showed that it could explain the solar and atmospheric 
neutrino deficits and give a dark matter candidate. We also discussed that there could 
be two typical parameter settings for Yukawa coupling constants, which brought about 
rather different phenomenological features. An interesing aspect of this model is that one 
of these parameter settings can also realize the mass and mixing pattern which has been 
known to explain the LSND results, simultaneously. It may be considered as another 
interesting feature of our mass matrix that it can explain both deficits of the solar and 
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atmospheric neutrinos without severely constraining the charged lepton mass matrix as 
long as it has no large off-diagonal elements. Although these features simply come from 
the extension of the parameter space due to the introduction of the new light sterile 
neutrino species, this kind of investigation can be considered to have a sufficient meaning 
to show the way for the extension of the neutrino sector. 

The introduction of the light sterile neutrinos is usually considered to be artificial. 
However, their appearence seems to be not so unnatural as shown in this paper if we 
take account of the generation structure of quarks and charged leptons and also assume 
the constrained nonrenormalizable superpotential. One of such simple and promising 
candidates may be the extra U(l) models coming from the Eq models inspired by the 
perturbative superstring, in which the group theoretical constraints on the Yukawa cou- 
plings are very weak. In such a case all but one or two of the right-handed neutrinos 
can be generally very light by the cooperation of both the superpotential constrained by 
the discrete symmetry and the extra U(l) D-fiat direction. They can play the important 
role in the neutrino physics such as neutrino oscillations. Although this scheme seems 
to be successful, it is generally not so easy to yield small neutrino masses and to induce 
the neutrino oscillations without bringing other phenomenological difficulties like proton 



decay in this framework pO|, 0. The simultaneous explanation of them will be the next 
step to build the realistic models in this direction. Anyway we believe that it will be 
worthy to proceed the further investigation of this kind of possibilities. 
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Table 1 Oscillation parameters for ■^j ijjfr. The positive contribution to the oscillation 
probability does not come from other combinations of a and (3. In this table we assume 
_ 2 i^j^g state identification V'/ = (-^a(= si), i^/^, Nb{= S2)) is assumed. 
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Table 2 Charge assingments under Zg x Zg discrete symmetry for lepton and Higgs 
sectors. Lf^ and Ej^ are the SU(2)2, doublet and singlet lepton chiral superfields, re- 
spectively. H12 and $ are the usual doublet Higgs and triplet Higgs chiral superfields. 
6 [i — 1,2,3) are the integers which can give the nontrivial charges to each fields and 
satisfy 1 < 6 ^ 8- 
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Table 3 Numerical values of neutrino oscillation parameters for the charged lepton 
mass matrices (A) and (B). We take parameters as a ~ 1.3, 6 ~ 2.4, cr ~ r ~ and 
^'e / ^¥ = 0.1. Assuming the mass differences which can explain the solar and atmospheric 
neutrino deficits, only the dominant contributions in each (■0/, ipf) sector are presented. 
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